Amorphous vortex glass phase in strongly disordered superconductors 
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We introduce a model describing vortices in strongly disordered three-dimensional superconduc- 
tors. The model focuses on the topological defects, i.e., dislocation lines, in an elastic description 
of the vortex lattice. The model is studied using Monte Carlo simulations, revealing a glass phase 
at low temperatures, separated by a continuous phase transition to the high temperature resistive 
vortex liquid phase. The critical exponents v w 1.3 and r\ ~ —0.4 characterizing the transition are 
obtained from finite size scaling. 
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Disorder has a profound impact on the phase diagram 
of type-II superconductors. The long range order of the 
Abrikosov vortex lattice in very clean samples is desta- 
bilized upon the introduction of disorder [l| and a glass 
state can appear instead The effects of weak dis- 

order are relatively well understood and may be treated 
within an elastic theory of the vortex lattice 0, 0, 0. 
The resulting phase, known as the Bragg glass or the 
elastic vortex glass, is characterized by an algebraic de- 
cay of translational correlations and perfect orientational 
order, and the absence of large dislocation loops Q. 

For stronger disorder dislocations start to become en- 
ergetically favorable and begin to proliferate. What hap- 
pens then is still a matter of controversy. One possibil- 
ity is that any kind of order in the system is lost and 
the system becomes thermodynamically equivalent to a 
vortex line liquid (although perhaps a very viscous one) 
with a non-zero resistivity [3j- The other possibility is 
that the dislocations themselves remain pinned and the 
system freezes into a different genuine thermodynamic 
amorphous glass phase, full of dislocation loops, and with 
zero linear resistivity. Such a disorder driven transition 
from the Bragg glass to the vortex glass have been pro- 
posed to occur as the magnetic field is increased 0, @ ■ 

Unfortunately, the theoretical methods based on pin- 
ning of elastic manifolds in random media, that were so 
successful for the Bragg glass, does not easily extend to 
incorporate the effects of dislocations. In this Letter we 
introduce and study a model that goes beyond a purely 
elastic description by including dislocations, and is well 
suited for numerical simulations. 

Experimentally, the evidence for a vortex glass comes 
mainly from the scaling relations obeyed by the current- 
voltage characteristics and resistivity at a continuous 
phase transition The first order vortex melting line, 
observed in very clean samples at low magnetic fields, is 
replaced by a smooth continuous transition upon increas- 
ing the field or disorder. The universal critical exponents 
extracted from measurements in this regime provide an 
important characterization of the glass transition that 
can be compared with theoretical predictions. A large 
number of studies have thus been able to obtain good 



scaling behavior in different materials, and this fact has 
been taken to support the existence of an amorphous 
glass phase 8J. This interpretation has, however, been 
questioned recently 0. In practice it may be extremely 
difficult to discriminate between a true thermodynamic 
glass and a very slowly moving non-equilibrium viscous 
vortex liquid (as window glass), with a small, but non- 
zero, resistivity due to thermally activated flux creep. 

On the theoretical side a lot of interest has focused 
on the so called gauge glass model This approach, 
which is inspired by spin glasses in random magnets, 
combines the U(l) symmetry of a superconducting or- 
der parameter with quenched disorder and frustration, 
which are believed necessary to create a glass phase, 
but has no other connection with the microscopies. The 
model, which consists of an XF-model in a quenched 
random gauge field, has a glass phase at low tempera- 
tures in three dimensions |l(1 | (but probably not in two). 
Once screening of the vortex interaction, due to gauge 
field fluctuations, is accounted for, the glass phase disap- 
pears 0. However, the way disorder enters the model, 
via quenched random fluxes instead of a random core en- 
ergy coupling to the vortices, is not very realistic. More 
importantly, perhaps, the model lacks a length scale, the 
inter-vortex distance, set by the external magnetic field. 

Attempts to study more realistic models have come to 
contradictory conclusions. Simulations of frustrated 3D 
Ay-modcls with quenched disorder show evidence for a 
direct disorder driven transition from the Bragg glass to 
a vortex liquid at low disorder and fillings, and find no 
clear signs of any amorphous glass phase [l2j (See, how- 
ever, Ref . Il3j) . On the other hand recent simulations 0] 
of a random pinning vortex model for much stronger dis- 
order and very high field, corresponding to a half-filled 
system, do find a continuous transition, but with criti- 
cal exponents quite far from experiments. None of these 
models included screening. Including screening in a sim- 
ilar model appears to destroy the glass phase 0], j us * 
as in the gauge glass. Langevin simulations of vortex 
lines also found no glass phase 0] . Recently, Kierfeld et 
al. 0] invoked a Landau theory in terms of the disloca- 
tion density, assuming the dislocation lines were directed, 
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and proposed a phase diagram containing an amorphous 
glass phase. Some support for a glass phase has also been 
found using the boson analogy |17| . 

The approach taken in this Letter starts from an elastic 
theory of the vortex lattice, but modified to include dis- 
locations. The basic degrees of freedom in our model is 
the two-dimensional displacement vector and a tensor 
field B that accounts for the dislocations. Since the clean 
Abrikosov vortex lattice usually has hexagonal symme- 
try we define the model on a stacked triangular lattice 
with N = L 3 sites and unit lattice constant, and with 
the displacements living on the vertices and the field 
By living on the links connecting nearest neighbors i and 
j. The Hamiltonian reads 



where (■••) denotes a thermal average, and the 3x2 
stress tensor a is given by 



Unfortunately, both of these vanish upon averaging over 
the disorder realizations, [(c)] av = [C] av = 0, since 
[/lav becomes independent of for the fully disordered 
model, where the averaging includes all possible deforma- 
tions. A non-zero (er) is still indicative of a phase with 
rigidity, and hence we use 
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where (ij) denote all nearest neighbor pairs on the 



instead, where the second form, which involves two differ- 
ent replicas a and /? to avoid any bias in the expectation 
value, is the one used in the actual simulations. Alter- 



stackcd triangular lattice, = u, - u., + B y -, and natively, of ms can be calculated from Eq. J2J using that 
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The By are integer multiples of the elemen- 
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tary triangular lattice vectors ai = (1,0), a 2 = (|> 2 ) 
in the plane, i.e., By = mai + n&2, m,n G Z, and con- 
tain all the information about the dislocations. More pre- 
cisely the dislocation density is given by the lattice curl of 
B, implying that dislocations form lines with conserved 
Burgers vector as they should. A dual representation in 
terms of interacting dislocation lines can be obtained by 
integrating out u, but Eq. turns out to be more con- 
venient for our purposes. The coupling to disorder occurs 
via the random variable 25]. For simplicity we take 
the maximum possible disorder strength, with F,y inde- 
pendently and uniformly distributed in the planar unit 
cell of the triangular lattice. 

In the absence of dislocations all B^ are zero and the 
model reduces to an ordinary elastic model of the vortex 
lattice with compression modulus en = 3/2, tilt modulus 
C44 = 1, and shear modulus cqq — 3/2 |2jj. The general- 
ization to other values of these (bare) elastic constants is 
straightforward via the addition of a term proportional 
to Yltij) ( w ij ' e ij) 2 i e ij being the unit vector connect- 
ing nearest neighbor sites i and j] to Eq. and/or 
by changing the coupling in the z-direction to make it 
anisotropic. In this work we will, however, restrict to the 
isotropic model defined above. 

The presence of thermally excited dislocations will 
renormalize the elastic constants and eventually cause a 
transition at which they drop to zero. The renormalized 
elastic constants can be obtained by taking the second 
derivative of the free energy with respect to a deforma- 
tion, Wjj — > w 
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In the solid low temperature phase (if it exists) ofms wn l 
acquire a non-zero expectation value, while in the liquid 
it approaches zero. Close to a continuous phase transition 
it obeys the finite size scaling form, 



Crms = L 1 d a(tL 



(6) 



where t = (T — T c )/T c is the reduced temperature, v the 
correlation length exponent, d — 3 the dimensionality, 
and a is a scaling function |27j . 

Lattice order may be probed by the structure function 
5k = jf (\^2i e%K ' Ui \ 2 \ at a reciprocal lattice vector K 
of the triangular lattice, but for strong disorder, when 
the system might freeze into a highly disordered glassy 
state, this is not useful. To quantify order in a glass one 
may instead consider an overlap order parameter, 
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between two replicas a and [3. The order parameter sus- 
ceptibility, given by 
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obeys the finite size scaling relation 
Xk = L 2 -"*(iLV") 
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which defines the critical exponent rj. 

We study this model via Monte Carlo (MC) simula- 
tions. One MC sweep consists of N randomly chosen 
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FIG. 2: Test of equilibration in the glass. As seen, af ms from 
FIG. 1: Clean system. The solid lines are interpolations Eq. @ and ® converge to a common value from below and 
using histogram reweighting j^. above, which provides a very useful check on the equilibration. 



heat bath updates of the and 12N Metropolis updates 
of By, where N = L 3 is the number of sites of the lattice. 
In order to overcome the slow glassy dynamics in the fully 
disordered system an exchange (or parallel tempering) 
MC algorithm is used, where several systems with 
identical disorder configurations but different tempera- 
tures are simulated in parallel and sometimes exchanged. 
The temperatures were chosen in a range around T c , such 
that the acceptance rate of temperature exchanges were 
no less than 15%. 

We begin with a brief discussion of the clean model 
(all Fij — 0). Figure n shows a sharp drop in the struc- 
ture function Sk and in the renormalized shear modulus 
C66, and a discontinuous jump in the average energy, with 
a corresponding approximate delta-function peak in the 
specific heat Cy . The very sharp variation in these quan- 
tities at T c k 0.076 signifies a discontinuous first order 
transition in accordance with experiments [l9j as well 
as numerical simulations of frustrated XF-models [2(| 
and vortex line models 0, • This is also confirmed 
by looking at histograms for the energy or the structure 
function, which show an increasingly growing double- 
peak structure developing for the largest system sizes. 

We now turn to the disordered model. In the maxi- 
mally disordered model the dynamics becomes very slow 
at low temperatures, and it is important to check the 
convergence of the MC simulation. To do so we mon- 
itor several quantities averaged over increasingly longer 
simulation times, after discarding an equally long time 
for equilibration. In Fig. [21 we plot, as an example, of ms 
calculated from Eq. J3J and JSJ), and xk as a function of 
MC time for the size L — 5 and the lowest T used in the 
exchange MC for this size. As seen the curves saturate 
and become time independent which indicates that the 
simulation has converged at around 100000 sweeps. Sim- 
ilar checks are done for all sizes. Typically, 4000-7000 
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FIG. 3: Finite size scaling collapse of MC data for ofms 
according to Eq. (JSJ. Corrections to scaling are clearly visible 
for the smallest size, L = 3, which was therefore omitted from 
the fitting procedure. The crossing point in the inset gives an 
estimate of T c . 



disorder realizations are used to form disorder averages. 

In the inset of Fig. |3] we show MC data for L 4 of ms 
vs. T for several different system sizes, which according 
to Eq. 10 should become independent of L exactly at 
T = T c . Indeed, the curves cross, which indicates that 
the model has a glass phase below this point, with a non- 
zero T c w 0.018. The main part shows a finite size scaling 
collapse of the same data, which gives an estimate of 
the critical exponent v. The structure function 5k stays 
constant 1 independent of temperature and system 
size, implying a total lack of crystalline order throughout 
the glass phase. Also, no signature of the transition is 
found in the density of dislocations, which remains finite 
down to the lowest temperatures. In Fig. 0] we plot a 
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study to include dynamic quantities, and to map out the 
phase diagram for intermediate disorder strengths and 
other values of the bare elastic constants. 
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FIG. 4: Finite size scaling collapse of MC data for \k using 
Eq. ©. 



finite size scaling collapse using Eq. © of xk (averaged 
over the smallest reciprocal lattice vectors) . As a measure 
of the quality of the data collapses we use the error in a 
X 2 -fit of the scaled data to a polynomial (indicated by 
the dotted lines in the figures). The critical parameters 
are estimated by adjusting them to minimize this error, 
and the errorbars are estimated by varying them until 
the data no longer scales. The final results for the critical 
parameters are summarized in Table [I] 

In summary, we have introduced a model of vortex 
lattices subjected to disorder. The model goes beyond 
an elastic description by including dislocations. Such 
a description should be valid as soon as the screening 
length is larger than the inter- vortex distance, A > a v . 
The model is studied via MC simulations and finite size 
scaling. Without disorder the model shows a first or- 
der melting transition. For strong disorder, the model 
has a low-temperature glass phase where the dislocation 
loops are pinned by the disorder. The glass melts via 
a continuous transition characterized by the critical ex- 
ponents in Table [I] We propose that the model gives 
an effective description of the vortex glass transition in 
type-II superconductors. Experiments usually find a cor- 
relation length exponent in the range 1 < v < 2 
which compares well with the present value v — 1.3±0.2. 
Very recent simulations on disordered 3D W-models 
find evidence for a glass transition with v = 1.5 ± 0.12 
or v = 1.1 ± 0.2, 77 = -0.5 ± 0.1 0, consistent with 
the present study. It would be interesting to extend the 



TABLE I: Estimates of the critical properties of the model. 
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This corresponds to a random stress. Other ways of incor- 
porating disorder would also be interesting to consider. 
To avoid factors of \/3/2 extensive quantities are normal- 
ized using N instead of the volume Q. = N\f3/2. 
The absence of any nontrivial anomalous dimension in 
Eq. © is a result of the underlying periodicity in the 
model. At a fixed point with such periodicity the scaling 
dimension of u must vanish, and hence w and B scales 
as L _1 and a as L 1_d . 



